Chiral logarithms to five loops 
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We investigate two specific Green functions in the framework of chiral perturba- 
tion theory. We show that, using analyticity and unitarity, their leading logarithmic 
singularities can be evaluated in the chiral limit to any desired order in the chi- 
ral expansion, with a modest calculational cost. The claim is illustrated with an 
evaluation of the leading logarithm for the scalar two-point function to five-loop 
, £t ■ order. 
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1 Chiral logarithms in the chiral limit 



In this article, we discuss the structure of chiral logarithms in the effective 
low-energy theory of QCD, chiral perturbation theory [1,2]. To simplify the 
discussion, we consider the case of two flavours u and d. The relevant effective 
Lagrangian is explicitly known at next-to-next-to-leading order in the chiral 
expansion [2-4]. Its structure is 

£ = £ 2 + £ 4 + £ 6 + ... , (1) 

where the indices stand for the chiral order of the Lagrangian. We count £ n a 
quantity of order p n . The low-energy expansion amounts to an expansion in 
powers of p 2 . The terms of order p 4 and p e e.g. are generated by 

£ 2 x £ 2 , £ 4 : p 4 , 

(2) 

£ 2 x £ 2 x £ 2 , £ 2 x £4, £ 6 : p 6 , 
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where the symbol x denotes a loop integral, and where C 2 x C 2 x £2 de- 
notes a two-loop integral, with 3 vertices from C 2 , etc. We furthermore use 
the convention that the low-energy constants (LECs) in C n ,n > 4 are made 
dimensionless by multiplying them with appropriate powers of F~ 2 , where 
F denotes the pion decay constant in the chiral limit. As usual, we call the 
terms of order p 2n+2 generically n-loop contributions, although they are not 
exclusively generated by n independent loop integrals. 

In the following, we consider the two-point function of two scalar quark cur- 
rents in the chiral limit m u = m d = 0, 

H(s) = i j dxe ipx (0\TS°{x)S°(0)\0) ; S° = uu + dd; s = p 2 , (3) 

whose low-energy expansion is of the form 

H(s) = H 2 + H 4 + H 6 + . . . , (4) 

where H n is considered to be of order p n according to the above counting. 
There is no tree graph contribution in this case, H 2 = 0. The leading term 
is [2] 

H 4 = B 2 [A 1 L + A ]; L = \n~ , 

A 1 = -JL,A = ± + 8(eZ + h{). (5) 

The quantity B is related to the quark condensate, (0|S' (0)|0) = —2BF 2 , and 
£3, h\ are renormalized LECs from £ 4 . The symbol L denotes a chiral logarithm 
which is generated by the relevant one-loop graph, and /i stands for the scale 
which is introduced using <i-dimensional regularization of the loop integral. 
The scale dependence of £3, h\ cancel the scale dependence of the chiral loga- 
rithm - the quantities H n are scale independent. 

Higher orders in the low-energy expansion generate additional powers of L. 
The two-loop contribution e.g. contains a term proportional to L 2 , 

sB 2 

H 6 = —[B 2 L 2 + BiL + B ] , (6) 

where 

B2 = 128^' (7) 
see below. In addition to the square L 2 , the two-loop term H 6 contains as well 
a single logarithm, as the one-loop contribution if 4 does. However, this term 
is suppressed here by one power in s, and becomes negligibly small at small 
momenta with respect to the one in H 4 . Likewise, the three-loop contribution 
H$ contains L n ,n < 3, where the term proportional to L 2 is again suppressed 
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Fig. 1. The power k of chiral logarithms of the correlator as a function of the 
number n of loops. The solid line corresponds to leading logarithms and the dashed 
line connects the terms of order I? . 

by one power of s with respect to the one in H 6 , etc. The general structure 
of the loop expansion is illustrated in Fig. 1, where we display the power of 
chiral logarithms as a function of the number of loops. 

In the following, we call the logarithm L n generated by n-loop graphs the 
leading logarithms (LL). They come with power s( ra_1 ) in the chiral expansion 
of H and correspond to the contributions connected by the tilted, solid line 
in Fig. 1. Terms of order L n generated by graphs with more than n loops are 
suppressed by additional powers of s, see the horizontal, dashed line in the 
figure, which connects the terms of order L 2 . 

Performing a renormalization group analysis, it can be shown [5] that the 
leading logarithm L n is generated by n-loop integrals of the type £ 2 x £ 2 x 
. . . x £ 2 (n + 1 vertices from C 2 ). Its coefficient does not, therefore, contain 
any LECs - it is determined by C 2 alone. 

There are two obvious questions: Is it possible 

i) to calculate the leading logarithm for any n? 

ii) to sum up the leading logarithms, similarly to summing up leading loga- 
rithmic singularities in renormalizable theories? 

Here, we concentrate on the first question, the second question will be ad- 
dressed in a subsequent publication [6]. There are several methods to evaluate 
the leading logarithms: 

i) Calculation by brute force. This procedure can obviously not be imple- 
mented beyond the first few terms. 

ii) The renormalization group can be used to show that the leading logarithm 
L n can be determined by a one-loop calculation with the Lagrangian L 2n 
for any n > 2 [5]. While this sounds promising, the method cannot, in 
practice, be used to perform the calculation beyond the three-loop level. 



3 



iii) The third method makes use of the fact that the correlator H(s) has 
a rather simple structure in the chiral limit. Unitarity, analyticity and 
the Roy equations then allow one to to determine the leading logarithms 
rather easily. 

We now illustrate methods ii) and iii) in turn, and start with ii). 



2 Leading logarithms from a one— loop calculation 



We illustrate the method with H 8 . We make use of the fact that H n are scale 
independent, 

H n = 0, (8) 
where we define the dot-operation by X = (J>^- After renormalization, in 
four dimensions, H 8 has the structure 



H* = B< 



+ C t L 



+ C 



+ d 2 (i;)l* + Di(£-)l + £>„(€/) 

+ Fi(cJ)L + F «) 
+ G (dl 



£ 2 x £ 2 x £ 2 x C 2 
C 2 x £ 2 x £ 4 

£4 x £4 

C 2 x £ 6 



(9) 

Here we have split the various contributions into different components which 
show the explicit dependence on the scale dependent LECs. According to 
the power counting rules it is possible to assign these components to the 
corresponding type of diagram. The polynomials -EVim-O^fc) fo r example stem 
from one loop diagrams with two vertices from £4. The c\ denote renormalized 
LECs from £ 6 . Scale independence demands 



c 3 = -mil 



mi) = \ [miiD + f 1 (4) + dm) - 4^ 



(10) 



In fact the last two terms vanish. C3 is determined by Ei(£-£l) and i*i(c[) 
which are given by one loop calculations with either two vertices from £ 4 
or one vertex from £ 6 and one vertex from £ 2 . Performing this calculation 
explicitly, one has to calculate the diagrams indicated in Fig. 2. The diagram 
Gc 4 xc 4 yields the contribution 



Gc^cAs) = -GtB 2 ^£jj(s) 



(11) 
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Gc 4 xc 4 Gc 6 xc 2 

Fig. 2. Diagrams used in the calculation of the three-loop leading logarithm of the 
scalar two-point function. The cross denotes a vertex from £2, the square denotes 
a vertex from £4 and the black dot stands for a vertex from £q. 

and for Gc 6 xc 2 one obtains 

Gc 6 xc 2 (s) = 9QiB 2 ^c 6 J(s), (12) 



with 

J(s 



16tt 2 



(13) 



4^ + r 'W- 1 H-4^J +2 J- 

Considering the coefficient of the logarithm L and using Eq. (10) one finds for 
the leading logarithm to three loops 

c --T(Is5^ + < s3 >- < 14 > 

Remark: The above argument shows that, to calculate C3, it is sufficient to 
know the Lagrangian £ 6 , whereas £g would be needed to renormalize all three- 
loop diagrams. Indeed, using the same arguments as above, it can be shown 
that a tree-level calculation with £§ suffices to calculate C3. Likewise, a two- 
loop calculation with £4 suffices as well to determine C3. End of remark . 

Whereas this method to calculate the leading logarithms is simple in princi- 
ple, it soon becomes prohibitively complicated: To evaluate the term L 4 by a 
one-loop calculation, one would need the renormalized £$■ We now present a 
method that allows one to calculate the LLs to arbitrary orders. 



3 Leading logarithms from unitarity and analyticity 



Here, it is useful to slightly rearrange the loop expansion of the correlator. We 
collect terms with the same power of the chiral logarithm and write 

H ( s ) = TT1 { P o + P i L + P2L 2 + ■■■}, (15) 

l07T Z 1 > 

where denote (dimensionless) polynomials in the variable iV = s/(167r 2 F 2 ), 
and in the scale dependent LECs. n-loop graphs contribute to polynomials 
P m with index m < n, with a term proportional to s^ n_1 \ Up to two loops the 
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leading contributions are given by 



P = 0, Pi = -6, P 2 = 6N. (16) 

In the following, we make use of the fact that H(s) is analytic in the complex 
s-plane, cut along the real positive axis. Unitarity of the S-matrix determines 
the discontinuity of H across this cut, 



H(s + ie) - H(s - ie) = ^(2vr) 4 ^ 4 )(P„, - j9)|(0|S°(0)|n)| 2 , 

n 

= ^\F(s)\ 2 + ..., (17) 
where F(s) is the scalar form factor, 

(0|S°(0)|7T^7rV)> = S lk F(s) , a = (p + pf. (18) 

The ellipsis in (17) denotes intermediate states with more than two pions. The 
main point is the observation that knowledge of the discontinuity allows one 
to reconstruct the leading logarithmic term easily, by constructing a function 
with the given discontinuity. This is trivial for a function with the structure 
displayed in (15). Assuming that only two pion intermediate states contribute 
to the leading chiral logarithrrJjO, the relation (17) allows us to calculate the 
TV-loop leading logarithm of the scalar two-point function, once the N—l loop 
leading logarithm of the scalar form factor is given. As the leading logarithm 
to two loops of the scalar form factor is known [7,8], the leading contribution 
of the polynomial P3 can be calculated along the lines just mentioned, 

Pa = -yiV 2 - (19) 

This result is in agreement with the relation (14) obtained from renormaliza- 
tion group arguments, and justifies the neglection of four pion intermediate 
states in the unitarity relation (17). In appendix B, a more general argument 
is given. 

This method allows one to calculate quite easily also higher-order leading 
logarithms. Indeed, as just seen, the scalar form factor F(s) determines the 
discontinuity of H(s). Unitarity may again be used to determine chiral loga- 
rithms in F(s), and hence in H(s). Let us illustrate the procedure by deter- 
mining the leading term in the polynomial P4 in (15), which amounts to the 
determination of the LL at four-loop order. The scalar form factor is analytic 
in the complex s-plane, cut along the positive real axis. Its chiral expansion 
can be arranged in the same way as for the correlator in Eq. (15), 

F(s) = 2B {t + T\L + T 2 L 2 + •••}, (20) 



This assumption will be justified below. 
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where the coefficients T« are again dimensionless polynomials in N and in the 
LECs. The leading contributions to the tree-level, one- and two-loop result 
are given by 

43 

T = l, T 1 = -N, t 2 = —N 2 . (21) 

The discontinuity across the cut of F(s) is given by 

F(s + ie) - F(s - ie) = 2i t° (s)F*(s) + .... (22) 

The quantity t° is the isospin 1 = S-wave of the 7T7r-scattering amplitude, 
in the notation of Ref. [2]. The ellipsis in (22) denotes contributions from 
intermediate states with more than two pions. Again discarding these for the 
moment, the relation (22) shows that knowledge of the LLs of both, the scalar 
form factor and t° at N loops, determines the LL of the scalar form factor at 
N + 1 loops, which then determines the LL in H(s) to iV + 2 loops. 

Here, we use the fact that t° is known up to two loops [9]. Together with 
the two-loop leading logarithm of F(s) [7,8], one finally finds the three-loop 
result for F(s) and the four-loop result for H(s), respectively, 

T 3 = -^iV 3 + 0(A P 4 = fN 3 + 0(s% (23) 

Again, we explicitly checked that T3 is correct by using the renormalization 
group equation of [5] as presented in Sect. 2. 



4 Leading logarithms to five— loop order 



The arguments in the previous section show that the correlator H, the scalar 
form factor F and the isospin zero S-wave determine a closed system as far as 
LLs are concerned. The chain is 

f° (s) ^ F(s) ^ H(s) , (24) 

where each step increases the knowledge of the LL by one-loop. The analytic 
structure of t° Q differs from the one of H and F. While the correlator H and 
the scalar form factor F are holomorphic in the complex s-plane cut along 
the positive real axis, this is not the case for t° . Indeed, this partial wave is 
holomorphic in the complex s-plane cut along the real axis in the intervals 
—00, 0) and (AMI, 00) for 

non vanishing quark masses. In the chiral limit 
considered here, the branch points at s = and s = coalesce: the partial 
wave contains, near s = 0, singularities of the type L_ = In— as well as of 
the type L + — In The latter are generated by the left hand cut. 
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As a result of this structure, unitarity alone does not provide sufficient infor- 
mation to determine the LL in t° . We illustrate the problem at the two-loop 
level: The structure of LL-terms of t° reads at this order 

t° (s) = jf{aLl + bL + L_ + cL\} + ■■■. (25) 

To calculate the discontinuity of F(s) according to Eq. (22), only the sum of 
the coefficients of all leading logarithms is required. It is determined by the dif- 
ference between the imaginary part for s = —\s\ + ie, Im~tg = j^iii(b + 2c)L + 

and the imaginary part for s — \s\ + ie, Ihi + £q = — J^7r(2a + b)L + . 
Using unitarity restricted to two pion intermediate states for t° 

+ .-, (26) 

only determines the imaginary part for s > 0. However, there exists a set 
of integral equations for the partial wave amplitudes of elastic tttt scattering, 
the Roy equations [10]. These equations provide the necessary information to 
determine the left hand singularity. 

This allows us to calculate the leading part of the polynomial P 5 . Let us switch 
on the quark masses for the moment. The Roy equations are of the form 

2 °° poo 

t\{s) = k\{s) +EE l 2 ds'K" (s, s')lm4 (s>). (27) 

J'=0f=0' ,4M ? 

where fcj(s) * s a subtraction polynomial and Kgf(s,s') are known functions. 
To determine the left hand cut, knowledge of the imaginary part of t° (s) 
for negative values of s is sufficient. One can therefore drop the subtraction 
polynomial k° (s) and one only needs the imaginary parts of the integration 
kernels as given in appendix A, which only differ from zero in a finite interval 
of s'. As the integration over s' only runs over positive values of s', unitarity 
provides the imaginary parts of the partial waves in the integrand, 



t° (s + ie) - t° (s - ie) = 2i t° (s) 



2 oo piMl-s 

Im-i°(s) = E E / 2 ds'lmKll, (s, s')lm + 4(s'). (28) 

As we are only interested in terms which generate leading logarithms we focus 
on parts of Imt^ which are proportional to L 2 + . 

We extend the function Imif °{,'(s, s') analytically to the whole complex s- 
plane and choose s = —AMI ~ $ + ^ e ^ > and e > 0. This allows to 
perform the chiral limit and still remain in a region where the real part of s is 
negative. In the infinite sum, only terms with £' < 2 contribute to the leading 
logarithm. Following once again the established path and using the unitarity 
relations several times, one obtains the four-loop LL of the scalar form factor 
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Fig. 3. Contribution of the LLs to the normalized scalar form factor in a region 
where F(s) is real. The solid, dotted, dashed and dashed-dotted lines denote the 
one-loop, two-loop, three-loop and four-loop LLs, respectively. The solid grey line 
indicates the tree-level result. 

T 4 and the five-loop LL of the scalar two-point function P 5 , 



Ta 



15283 
9720 



N A + 0(s 5 ), 



140347 
16200 



iV 4 + 0(s 5 ). (29) 



We checked the sum of the coefficients of the leading logarithms of t° to the 
order p 8 with the renormalization group [5] . 

We shortly mention a technical complication that arises while using the Roy 
equations in the chiral limit. According to Eq. (28), the integrand is propor- 
tional the product of the imaginary parts of the Roy-kernels and of the partial 
waves. The centrifugal barrier disappears in the chiral limit, as a result of 
which the imaginary parts are all of the same order at threshold, independent 
of the angular momentum involved, 



limbs') 



~ s 



i 4 



i'>2 



On the other hand, the imaginary parts of the kernels behave like 



ImX//'( S , S ')~-(4 



(30) 



(31) 



' > 5, the integrand thus becomes infrared singular 
' = 5, a divergence of the form s A \n(M n ) is generated. 



in the chiral limit. For 
and non integrable. For 
Divergences of the form s e _1 M~ 2 ^ also occur - these are, however, at least 
of order 0(s 5 ) and do not affect the four-loop leading logarithm T 4 . The 
divergences never appear in the terms relevant for the calculation of the leading 
logarithms and do not, therefore, affect their coefficients. 

To illustrate the size of the LLs, we plot in Fig. 3 the numerical contributions 
of the LLs of the scalar form factor for two choices of the scale ji. One sees that 
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for s > —0.85 GeV 2 , the LL corrections of higher-orders are always suppressed 
with respect to the lower orders. 



5 Pade approximants 

Pade approximants may be used to estimate higher-order terms. The Pade 
approximant [M, N] for the scalar form factor is 

[M, N] = ^ = |M + (p^ N ^), (32) 

where Ek only contains terms up to and including order p 2K . In Ref. [7], 
it is shown that the Pade approximant [0, 1] does not reproduce the correct 
factor of the two-loop LL. This shortcoming persists for the factors of the 
higher LL. Also the Pade approximants of the two- and three-loop LLs, [1, 1], 
[0,2], [2,1], [1,2] and [0,3], fail to reproduce the correct three- and four- 
loop LLs, respectively. In Tab. 1, we compare the two- and three-loop Pade 
approximants with the exact result. 

Table 1 

Pade approximants for the scalar form factor. Displayed are the exact and Pade 
values for the coefficients of the leading logarithms. 





F(s) 
F(0) 


[0,1] 


[1,1] 


[0,2] 


[2, 1] 


[1,2] 


[0,3] 


TV 3 

T 4 
TV 4 


143 
108 

15283 
9720 


-1 
1 


1849 
1296 

79507 
46656 


25 
18 

2101 
1296 


143 
108 

20449 
13932 


143 
108 

1961 
1296 


143 
108 

1933 
1296 



6 Summary and conclusion 

1. In this letter, we study the structure of chiral logarithms in the two-flavour 
sector, in the chiral limit m u = = 0. In particular, we investigate the chiral 
expansions of the scaler form factor of the pion, and of the correlator of two 
isoscalar quark currents. We confine our interest to the so called leading chiral 
logarithms, which are defined to be the ones accompanied with the least power 
of external momenta. 

2. Standard techniques to calculate these logarithms - direct evaluation by 
calculating Feynman diagrams, or using renormalization group techniques - 
cannot be applied beyond three-loop accuracy, because the required labor 
becomes prohibitive. 
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3. We point out that, using unitarity and analyticity, the leading logarithms 
for the mentioned correlators can be calculated to any desired order with a 
modest calculational cost. We illustrate this claim with an evaluation of the 
leading logarithm for the scalar two-point function to five-loop order. As far 
as we are aware, this is the first evaluation of chiral logarithms at this accuracy. 
The numerical structure of the coefficients appears erratic to us. 

4. The proposed technique makes use of Roy-equations for the titt amplitude 
near threshold, in the chiral limit. While the Roy-kernels become singular in 
this limit, we argue (and explicitly check in one case) that these singularities 
do not affect the evaluation of the leading logarithms. 

5. We compare our results with Pade approximants and confirm an earlier 
statement [7] concerning the use of this technique for the evaluation of chiral 
logarithms: knowledge of the coefficients up to and including those of order N 
does not allow one to calculate the ones at order N + 1 and higher. 
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A Roy equations 

Using the integral representation for the kernels K^,'(s,s') given in [11] it is 
straightforward to calculate their imaginary parts, 

2s + s' -4M*\ 
4M2 - S ' J ' 

(A.l) 

where are the matrix elements of the crossing matrix C su and P n (x) are 
the Legendre polynomials. Note that it is crucial to use the projection interval 
[0, 1] in the integral representation of the kernels as given injll] . Otherwise, 
the Roy equations cease to be valid for negative values of so 



4 We are indebted to H. Leutwyler for discussions on this issue. 



t kW( * r iv 4f + 2 p ( s_ + 2s^-ml \ ( 
Im K u , (,,,) = C su j^Pi ( AM ,_ s 1 P* [ 
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B Suppression of higher intermediate states 

In this appendix, we argue that intermediate states with more than two pions 
do not contribute to the leading logarithms of the scalar two-point function. 
In the sum 

H(s + ie) - H(s - ie) = iJ2(2K)^ ) (Pn-p)\(0\S°(0)\n}\ 2 (B.l) 

n 

appear the matrix elements of the 2k pion intermediate states 

Y 2kw = (0\S°(0)\2kn). (B.2) 

The tree-level results of these matrix elements start with i^~ 2 ( fe-1 ). Because the 
integrand of the phase space integral is the modulus squared of Y 2 kn, the lowest 
order contribution of the 2k pion intermediate states to the discontinuity of 
H(s) is therefore of order F -4 ^ -1 ^ and does not contain any logarithms. The 
loop corrections to these matrix elements are suppressed even stronger in 
F~ 2 . Take for example k = 2. The four particle phase space integration would 
have to generate a logarithm squared to contribute to the leading term in the 
polynomial P3. This is not possible and one way to see it is the following: Only 
diagrams with vertices solely from C 2 contribute to the leading logarithms. 
The highest power of the pole term e _1 is directly related to the highest power 
of the logarithm. Since the phase space integration in an infrared save theory 
does not produce a divergence, a phase space integration does not increase the 
power of the pole term and thus neither the power of a logarithm. In the case 
of the scalar form factor, the argument is the same. The higher intermediate 
states are also suppressed in F~ 2 . 
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